arXiv: 1503.07169vl [hep-ph] 24 Mar 2015 


YITP-SB-14-23 


Strong Optimized Conservative Fermi-hAT Constraints on Dark 
Matter Models from the Inclusive Photon Spectrum 

Andrea Massari/’[^ Eder Izaguirre,^’[^ Rouven Essig,^’[^ 

Andrea Albert,^ Elliott Bloom,^ and German Arturo Gomez-Vargas^’® 

^C.N. Yang Institute for Theoretical Physics, 

Stony Brook University, Stony Brook, NY 11794-3840 
'^Perimeter Institute for Theoretical Physics, Waterloo, ON N2L 6B9, Canada 
^SLAC - KIPAC, Stanford University, Stanford, CA 94305-4085 
^Instituto de Fisica, Pontificia Universidad Catolica de Chile, 

Avenida Vicuna Mackenna 4860, Santiago, Chile 
^Istituto Nazionale di Fisica Nucleare, 

Sezione di Roma “Tor Vergata”, 1-00133 Roma, Italy 

Abstract 

We set conservative, robust constraints on the annihilation and decay of dark matter into various 
Standard Model hnal states under various assumptions about the distribution of the dark matter in 
the Milky Way halo. We use the inclusive photon spectrum observed by the Fermi Gamma-ray Space 
Telescope through its main instrument, the Large-Area Telescope (LAT). We use simulated data to 
hrst hnd the “optimal” regions of interest in the 7-ray sky, where the expected dark matter signal is 
largest compared with the expected astrophysical foregrounds. We then require the predicted dark 
matter signal to be less than the observed photon counts in the a priori optimal regions. This yields 
a very conservative constraint as we do not attempt to model or subtract astrophysical foregrounds. 
The resulting limits are competitive with other existing limits, and, for some hnal states with cuspy 
dark-matter distributions in the Galactic Center region, disfavor the typical cross section required 
during freeze-out for a weakly interacting massive particle (WIMP) to obtain the observed relic 
abundance. 
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1. INTRODUCTION 


The Fermi Gamma-ray Space Telescope (Fermi), through its main instrument, the Large 
Area Telescope (LAT) |T], has been surveying the 7 -ray sky since August 2008 in the energy 
range from 20 MeV to above 300 GeV (with detected events up to ~ 1 TeV). In addition to 
7 rays produced by known astrophysical sources, the Fermi-LAT can detect photons from 
postulated decay or annihilation of dark matter (DM) to Standard Model (SM) particles. 
The possibility that DM can annihilate is particularly motivated by the “WIMP miracle” [2]. 
Here one hypothesizes the existence of weakly interacting massive particles (WIMPs) with 
few-GeV to few-TeV masses and weak-scale annihilation cross sections. These WIMPs would 
have been in thermal equilibrium with the SM sector in the early Universe and they generally 
produce the observed relic abundance of DM from thermal freeze-out. This suggests that 
WIMPs could still be annihilating today to SM particles. The annihilation could produce 
various SM particles, which can either radiate photons, further decay to other SM particles 
including photons, or inverse Gompton scatter (IGS) off background light, producing high- 
energy 7 rays. Those photons that arrive at the Fermz-LAT could then be used to infer 
properties of the DM particles and their distribution around us. 

Many WIMP searches have been performed using Fermz-LAT data. Analyses by the 
Fermi-LAT Gollaboration and outside groups have searched for monochromatic 7 -ray 
lines |3H9] and continuum 7 -ray excesses in the diffuse spectrum from different target 
regions e.g., dwarf spheroidal galaxies [ 101116 ] . clusters of galaxies [TTHI^ . the Galactic 
halo [131 120H22] . the Inner Galaxy |23ll36] . the Smith cloud jST] [38], and the extragalactic 
7 -ray background [39ll42j . No undisputed signal of DM has been detected thus far, and the 
cross-section upper limits from these analyses for DM masses uidm < 10 GeV are approach¬ 
ing the typical cross section required during freeze-out for a WIMP to obtain the observed 
relic abundance, namely {cru)reiic ~ 3 x 10“^® cm^ s“^. 

While DM is often thought of as being a stable particle, viable DM candidates only 
need to be stable on cosmological time-scales. In particular, DM lifetimes of the order 
of the age of the Universe or longer (tdm > 10^^ s) can typically evade cosmological and 
astrophysical bounds more easily than annihilating DM, such as constraints from Big Bang 
Nucleosynthesis |13], the extragalactic 7 -ray background m, and re-ionization and the 
Gosmic-Microwave-Background [45ll50j . The more relaxed constraints on decaying DM are 
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a result of the DM decay rate being linear with poM, as opposed to quadratic with pdm in 
the case of annihilation. 

In this paper, we will provide conservative DM cross-section upper limits and decay¬ 
lifetime lower limits from the Fermz-LAT inclusive photon spectrum. The inclusive spectrum 
is presumably dominated by astrophysical foregrounds in the Milky Way, though DM could 
contribute to it. We make no attempt at subtracting foregrounds and simply require that 
any putative DM signal contribute less than the observed flux. A similar idea has been 
used in other papers to derive conservative constraints 1201 Ea 123], where the DM signal is 
maximized until saturating the observed flux. The approach in this paper differs from such 
previous analyses in several ways, resulting in stronger constraints on DM. Firstly, we restrict 
our regions of interest (ROIs) to have a particular symmetric shape determined by only a 
few free parameters, and we optimize over these parameters. Secondly, we also optimize the 
energy range that we use for deriving the constraint. Thirdly, we optimize with respect to 
the constraint itself and not, for example, the signal-to-noise ratio, and last, we optimize our 
constraints on 10 simulated data sets, not on the measured data. After hnding the optimal 
ROI on simulated data, we use the real data from that same ROI to hnd the constraint. We 
derive constraints in this fashion for various DM-halo shapes and for various annihilation 
and decay hnal states. The resulting constraints, while being robust and conservative as 
no foregrounds have been subtracted, are competitive with other existing constraints and 
stronger than other conservative bounds obtained by 1201 [221 E3]. 

The paper is organized as follows. In ^ we discuss the calculation of the expected 7 - 
ray flux from DM annihilation and decay. In ^ we discuss the event selection, method, 
simulated data sets, and ROI selection. Q discusses the resulting constraints, while our 
conclusions are in ^ In Appendix we use our method to calculate the limits on DM- 
annihilation models that have been invoked to explain an excess of 7 rays from the Galactic 
Center (GC) and Inner Galaxy region. Appendix [^presents the optimal ROIs together with 
the corresponding count spectra for several DM channels. Appendix discusses the effect 
on our results of source masking and choosing front-/back-converting events. Appendix [Pj 
describes the astrophysical assumptions affecting the results that include contributions from 
ICS. Finally, Appendix [^provides more details on the simulated data sets that we use, and 
Appendix compares the limits obtained from our simulated data sets with those derived 
from real data. 
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2. EXPECTED DARK MATTER SIGNAL 


Gamma rays from DM annihilation or decay to SM final states can be produced in two 
dominant ways. The hrst possibility, which we refer to as prompt, is from either hnal-state 
radiation (FSR) produced by Bremsstrahlung by SM particles or from the decay of hadrons 
that arise in hadronic hnal states. The second possibility is from electrons and positrons 
(produced either directly or at the end of a cascade decay chain) that inverse Compton scatter 
off background ambient light, which primarily consists of starlight, the infrared background 
light, and the Cosmic Microwave Background (CMB). This ICS process boosts the energy 
of the background light to produce 7 rays. Unlike prompt radiation, ICS depends on various 
unknown astrophysical parameters discussed below. Although a sizable contribution to the 
energy lost by the electrons propagating through the Galaxy consists of synchrotron radiation 
due to acceleration by the Galactic magnetic held, we note that the synchrotron radiation 
does not make up a noticeable fraction of the 7 rays in the energy range under study, as 
we only consider DM particles with mass below 10 TeV [CT [52] . We thus do not include it 
in this study. Moreover, the DM signal can receive additional sizeable contributions due to 
Galactic substructure, particularly for annihilations [53], but we do not include this effect in 
our study. This makes our analysis more conservative and model independent in this regard. 
We now outline the calculation of the DM-initiated 7 -ray hux. 

2.1. Prompt radiation 

The differential hux, d<h.y/dU..y, of prompt photons coming from DM annihilation within 
the Milky Way halo is given by 

2 j 

where (crn) is the thermally averaged DM annihilation cross section, tudm is the DM mass, 
and (\.N^/dE^ is the photon spectrum per annihilation. We assume p© = 0.4 GeV/cm^ is the 
DM density at the Sun’s location in the Galaxy |5ll [55] and Tq = 8.5 kpc is the distance 

^ A range of values between 0.2 and 0.85 GeV/cm^ are possible at present though |54t[58] . Note that a 
different value for the local DM density would shift up or down our predictions for DM annihilation and 
decay by a factor proportional to and poj respectively for annihilations and decays. 
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between the Sun and the GC |59]. The “J-factor” is given by 
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which depends on the distribution of DM in the Milky Way halo, p{r), where r = r{s,b,i) 
is the Galactocentric distance, given by r = + Tq — 2srQ cos i cos 6, where i and b are 

the Galactic longitude and latitude, respectively, and s is the line-of-sight distance. The 
integral is over a particular ROI. For decays we can replace {av)p^/2m^y^ with pQ/rmoM 
in Eq. Q, where r is the DM decay lifetime, with the J-factor 

(3) 

JKOl ^0 P© 


Moreover, for decays the dN^/dE^ should be interpreted as the photon spectrum for indi¬ 
vidual DM particle decays. 

We consider four popular DM density profiles: the Navarro-Frenk-White (NFW) |6ni l6T] . 
Einasto |62l[63]. Isothermal [Mjj^and “contracted” NFW (NFWc) [MlET] with slope values 
taken from j5]. 
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We set a = 0.17, = 20kpc [^31 EZ], and = 5kpc jS3]- The normalization p{rQ) = Pq 

fixes Pq° ~ 1.56, Pg^^ — 0.35, pg™ ~ 0.08, and Pg^^"" ^ 0.24 in units of GeV/cm^. Our 

choice of the functional form and parameters in Eq. ([^ is a representative example of the 

possibility that, due to adiabatic contraction from the inclusion of baryonic matter, the DM 

profile might have a central slope steeper even than that of the NFW or Einasto profiles 

(although note that high-resolution observations of the rotation curves of dwarf and low- 

surface-brightness galaxies favor cored distributions 169]). The four profiles are shown 

^ Another popular parametrization of a “cored” profile is the Burkert profile jB^l. Adopting the best-fit 
parameters in [55] yields a distribution that is very close to the Isothermal one for radii < lOkpc. We will 
see that the optimal ROIs for the Isothermal profile for DM annihilation are contained with this region, 
and therefore the limits for the two cored distributions would be very similar. Thus we not explicitly 
consider the Burkert profile in our analysis. c 
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Figure 1. Left: Dark-matter density profiles versus distance from the Galactic Center (GC). We 
use the Isothermal (green), NFW (red), Einasto (blue), and a “contracted” NFW (NFWc, orange, 
with p oc 1/r^'^ for r —)• 0) prohle. Right: Prompt 7 -ray spectra produced in the annihilation of 
1 TeV dark matter to r+T“, bb, W^W~, uu, gg {g = a gluon), and 00, where 0 

decays either only to e+e" (with = 0.1 GeV), or only to p~ (with = 0.9 GeV), or to 
e'*“e“, and tt+tt" in the ratio 1:1:2 (with = 0.9 GeV). 


in Fig. (left). 

The (prompt) photon spectra, dN^/dE^ have been generated with Pythia 8.165 [70] 
or are based on formulas in [TltiTl] . They are the same as in DMFIT [75] after the latest 
update described in na. We will consider the ten different final states e’''e , , 

bb, W~^W~, uu, gg {g = a gluon), and 00 , where 0 decays either only to (with 

= 0.1 GeV), or only to p~^p~ (with = 0.9 GeV), or to e+e“, p~^p~, and tt+tt" in 
the ratio 1 : 1:2 (with = 0.9 GeV) (the latter ratio is motivated if 0 is a dark photon 
that kinetically mixes with the SM hypercharge gauge boson). Other SM hnal states are of 
course possible but they would yield constraints very similar to the channels we consider in 
our analysis. The annihilation channels to 00 are motivated by DM models [76l|77| that 
attempt to explain the rising positron fraction measured by PAMELA [78], Fermi |79j . and 
AMS-02 [Snil81j : the 0 can also facilitate an inelastic transition between the DM ground state 
and an excited state [zniisi to explain e.g., the 511 keV line anomaly [S3j. For DM decays, 
the 0 channels can be viewed as “simplihed models” that can capture how the constraints 
change when there is a cascade, e.g., [52]. We will sometimes refer to these scalar-mediated 
processes as “eXciting Dark Matter” (XDM). These spectra are shown in Fig. (right) in 
the case of annihilating DM and ttidm = 1 TeV. We do not consider other popular DM 
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candidates like axions and gravitinos. 


We note that the observed differential photon flnx can also be written as 


d$^ _ diVf), _ 1 diV^ 

ttot ^eflf dE^ £ dE^ ’ 

where we have now explicitly included Aes, the effective area (which is a function of energy), 
ttot, the LAT’s total live time, and £, the LAT’s exposure. Given the photon spectra, the 
number of photons from a DM annihilation signal in a spatial region D®, with J-factor 
and energy range [E^, is given by 
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where is the exposure averaged over D* and calculated at the midpoint of [Ek, Ek+i] 
(since the variation of the exposure over a single energy bin is very small). For decays the 
predicted counts are 



1 ^0P0 1 ji ci,k 
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The approximately homogeneously distributed DM in the Universe could provide an extra- 
galactic contribution to the observed photon flux. However, the observed y-ray spectrum 
will be different than that expected from Galactic DM interactions since the photons redshift 
as they propagate to us and there is a finite optical depth — the result of interactions of 
the 7 rays with low-energy photons that compose the extragalactic background light (EBL). 
This yields the following expected extragalactic photon intensity for decaying DM |M1185] 


d^d)., _ 1 DdmPcO r, diV 

dU^^dU dvr rmDM Jo H{z) dEy ^ 


Here, Ddm — 0.267 is the present DM energy density, Pc,o — d.7 x 10“® GeV/cm^ is 
the critical density today, Ey{z) = Ey{z + 1) is the energy of the emitted photon, 
H{z) = Ho\/^£lmA + + ^A, where Aim — 0.317 and Da — 0.683 are the total matter 

and cosmological-constant energy densities jHH], respectively, and we assume a flat Universe 
with Dm + DA = 1. The optical depth is given by r {Ey, z), and we use the parameterizations 
found in |8l]. We note that, for annihilating DM, the smooth extragalactic contribution is 













sub leading compared to the Galactic one and we ignore it, whereas for decays it is a factor 
of order < 1 as large as its Galactic counterpart and we include it in our analysis. 


2.2. Inverse Compton Scattering 


We include the flux generated by IGS for the cases where DM annihilates/decays to 
, as well as 00 channels. In all cases we end up with high-energy electrons 
and positrons. These propagate within the Galaxy and can lose energy through IGS off 
starlight, infrared background light, or GMB photons, or via synchrotron radiation in the 
Galactic magnetic held. The IGS process (e^'y' —)■ e^y) can produce high-energy y rays 
that are observed by the LAT. The synchrotron-radiation contribution in the Fermi-LAT 
energy range is subdominant for the DM masses and models that we consider [52] and thus 
we neglect it when we derive our limits. However, it must be included when calculating the 
IGS y-ray signal, since it determines how fast the electrons and positrons cool. In fact, the 
cooling time is strongly dependent on the Galactic magnetic held, whose values at diherent 
locations in our Galaxy are not known very accurately. This leads to large uncertainties in 
the IGS signal. Moreover, the calculation of the IGS signal requires additional assumptions; 
for example, we assume, as generally done, that the density of electrons and positrons 
after propagation follows a steady-state solution. However, phenomena such as the Fermi 
bubbles |87) . pointing to a dynamical event in the Milky Way’s recent history, might make 
this assumption not fully justihed. We also assume that the steady-state propagation of the 
electrons/positrons only occurs inside a cylindrical region of the Galaxy that has a maximum 
radius Rh and half-height Zh- The steady-state dihusion equation is given by (e.g., [58] ) 




dE' dE'^ 
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Here drie/dE'^ = dne{r,z,E'^/dE'^ is the energy-dependent diherential electron+positron 
density at a given point in the Galaxy, (r, z), where r and z are the cylindrical coordinates 


of the electron/positron in the Galaxy. The right-hand side of Eq. (12) is the source term 
and contains the DM density prohle, p{r, z) (a function of cylindrical coordinates) and 
the electron+positron energy spectrum, dW/dE'; also, there is a factor 1/2 for Majorana 
fermions, otherwise 1 for Dirac fermions. The hrst term on the left-hand side accounts for 
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the spatial diffusion and is characterized by an energy-dependent coefficient, 


(K) 



( 13 ) 


The second term is the energy-dependent loss and is given by 


b {r, z, E') = 


dEL 4a- 


dt 


3m2 




UB{r,z) + ir,z)Ri{E'^) 


2 = 1 


(14) 


where ar = Swr^/S, with re = aem/fne-, is the Thomson cross section, and ub (r, z) = B‘^/2 
is the energy density of the Galactic magnetic field B, chosen to have the form |89] 


B = B{r,z) = , 


(15) 


where Rb = 10 kpc and Zb = 2 kpc. The u-yj (r, z) are the energy densities of the three relevant 
light components in the Galaxy, i.e.: GMB, infrared light, and starlight. The factors Ri{E'^ 
take into account relativistic corrections. The y-ray differential flux at energy E^, resulting 
from IGS off an electron is 


di?y,dr2 


a: 


ds 
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where s is the line-of-sight distance, and 
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(18) 


We calculate the IGS contribution with GALPROP V50 [SO]. We use a version of GALPROP V50 
that was modified by the authors of |91) to include various DM annihilation and decay final 
states. We fix 5 = 0.33, Eq = d GeV, and take the cylindrical geometry to have a maximum 
radius Rh = 20 kpc and a maximum half-height Zh = d kpc. As mentioned above, the 
greatest source of uncertainty is due to the Galactic magnetic field, B. To capture some of 
this uncertainty, we vary Bq between 1 — 10 /iG, when showing our results in ^ We fix the 
spatial diffusion coefficient parameter to be Dq = 4.797 x 10^® cm^/s (6.311 x 10^® cm^/s) 
for Bo = 1 fxG (10 yuG). (See Appendix [Pj for sources for these values.) In Appendix [Pj we 


10 














show how our results are affected when varying Zh and Rh, in addition to Dq and Bq] we 
hnd that the largest effect on the results comes from the variation of Bq. 


3. DATA SETS AND METHODS 

We aim to set conservative, robust constraints on the annihilation and decay of DM into 
various SM final states. We consider the inclusive photon spectrum observed by the Fermi- 
LAT, and use simulated data to first find the “optimal” ROI in the 7 -ray sky, i.e. the one 
that yields the strongest constraint. We then require the DM signal to be less than the 
observed photon counts. We note that our approach does not allow us to search for the 
existence of a DM signal. 

In this section we describe the event selection, how we use the simulated data sets in 
our analysis, the ROI choice, and how we construct optimal upper bounds on the DM 
annihilation cross section and lower bound on the DM decay lifetime. We also provide a 
detailed example of our procedure. 


3.1. Event Selection 

The data set used for this study consists of ~ 5.84 years of Fermz-LAT data (from 
August 2008 until June 2014) in the energy range 1.5 — 750 GeV. We select photons using 
the P7REP_CLEAN event-class selection |92], to minimize contamination by residual cosmic 
rays. We also require the zenith angle to be smaller than 100° to remove photons originating 
from the bright Earth’s Limb. Details on the Fermi-LAT instrument and performance 
can be found in lU [93]. All data reduction and calculation of the exposure maps were 
performed using the Fermi-LAT ScienceTools, version v9r34pl |94| . As for the Fermi-LPFL 
instrument response functions (IRFs), we use P7REP_CLEAN_V15 for both MCs and data. 
As described in Appendix the results shown in this paper are obtained after masking all 
known point sources identified in the 5-year Fermi catalog (3FGL) [93], using a PSF (point 
spread function)-hke masking radius, except for those photons coming from within the inner 
2° X 2° square at the GG. Moreover, we include both front- and back-converting events. In 
Appendix we show that, although this choice is generally optimal, our results are not 
significantly affected if we mask only the brightest sources, or no sources at all, and if we 
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include only front- or only back-converting events. 


3.2. Simulated (Monte Carlo) Data Sets 


For our study, we use 10 Monte Carlo (MC) data sets, each a statistically independent 
~ 5.84-year representation of the 7 -ray sky. The same event selection described above is 
applied to MC data. We use the simulated data sets to select “optimized” ROIs, indepen¬ 
dent of the real data, as described below in §3.3 By hnding optimal ROIs based on the 
MC simulations, we avoid the possibility of accidentally obtaining a strong constraint due 
to statistical fluctuations in the data. We describe the details of the simulated data in Ap¬ 
pendix Note that the MC simulations contain photons with an energy range of 0.5 GeV 
to 500 GeV (as opposed to 1.5 GeV to 750 GeV in the data). We account for this difference 
by extrapolating the MG data up to 750 GeV as described in Appendix |E| 


3.3. ROI Choice 


We take the ROI for annihilating DM to have the dumbbell shape as shown in Fig. 
(left). This shape depends on three parameters: the radius from the GG to the edge of the 
ROI, R, the width in latitude of the Galactic Plane (GP) that is to be excluded from the 
ROI, 2Ab, and the width in longitude of the GG region that is to be included in the ROI, 
2Ai. The motivation for choosing such shape is that the DM distribution is approximately 
spherically symmetric (hence the choice of a circular region, parametrized by R), but the 
Galactic foregrounds are largest in the GP region, which we then remove. However, we 
include the GG in our ROI as this is where the DM signal peaks as well, dramatically so for 
cuspy prohles (since A^,dm oc Pdm)- 

For decaying DM, our choice of ROI will consist of the two high-latitude regions shown in 
Fig. 1^ (right), and depends on only one parameter: the width in latitude from the Galactic 
poles to the edge of the ROI, Ab^. In contrast to annihilation, the decaying DM signal is 
expected to be much less concentrated in the GG, since A^,dm cc Pdm- 
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Figure 2. Left: The choice of ROI (shaded) in the 7 -ray sky for dark-matter annihilation. The 
ROI depends on 3 parameters, as indicated. Right: The choice of ROI for dark-matter decays 
(shaded), which depends on one parameter, as indicated. 


3.4. Optimizing the ROIs and Energy Ranges using Simulated Data 


A particular DM model or Theory Hypothesis, Th = [utdm, P, annihilation/decay, channel], 
is characterized by the DM mass (mDivi), the DM density prohle (p), whether it is annihi¬ 
lating/decaying DM, and the annihilation/decay hnal state. Given any ROI and a photon 
energy range, [ROI, AE], we obtain a constraint on either the DM annihilation cross section, 
(an), or decay lifetime, r, for a given Th by requiring that the number of DM events, 
in [ROI, AE] does not exceed the observed value, N^^o- More precisely, to set a limit with 
a conhdence level (C.L.) of 1 — a, we vary (an) or r until the probability that > Af-y^o 

is a; in equations, the bound on (an) or r is obtained by solving 


^~f,0 

Poisson {k I ,dm) — 1 — a, (19) 

k=0 

where as usual 


Poisson(P I A) 


e ^ 
k\ 


( 20 ) 


For each Th, we hnd the optimal ROI and optimal photon energy range, [ROI, AE]o, which 
provides the best limit on (an) or r. If we simply scan over all [ROI, AE] in the data, this 
would subject our constraints to statistical fluctuations. Instead, we use the 10 simulated 
data sets to hnd [ROI, AE]o as follows. For the Tth ROI and energy range, [ROI, AE]i, 
and j-th simulated data set, we calculate the bound on the cross section or lifetime, (an)jj 


or Tij, as described above in Eq. (19). We then average the resulting expected limit across 
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the 10 simulations, i.e. 
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( 21 ) 

( 22 ) 


j = l 

We then hnd [ROI, AE]o by scanning over all [ROI, AE]i's and selecting the one that yields 
the minimum (uu)- (maximum Xj), i.e. 


{av) = min {(Jv )-, 

i 

T = max Ti . 

i 


(23) 

(24) 


We then use [ROI, AE]o on the real data to calculate the limits on (av) or r for the given 
Th- 


The ROIs used in our optimization are given in ^3.3 We bin each simulated data 
set into 0.18° x 0.18° rectangular pixels in Galactic latitude and longitude and N = 127 
logarithmically-uniform energy bins between 1.5 — 750 GeV. We then vary the ROI shape 


parameters described in ^3^ in steps of 0.5° for R, steps ~ 1° for Ab and Ai, and 2.5° for 
Abd- For each choice of ROI, we scan over all {N — 1){N — 2)/2 = 8064 choices of adjacent 
bins in energy, assuming a minimum of three adjacent bins. In Appendix we show a 
sample of the resulting optimized ROIs and energy ranges. 


We note that for large enough or the statistical distributions are approx¬ 
imately Gaussian, and we would obtain a 95% G.L. bound by requiring < Af^,o + 

1.64-y/iV^ ~ N^^o- Even our smallest optimal ROIs with the highest optimal energy ranges 
contain at least 0(10) photons. Our method thus does not produce constraints that are 
susceptible to Poisson fluctuations of the number of events in [ROI, AE]o, and, as a con¬ 
sequence, our constraints are not expected to improve signihcantly with more data (some 
small improvements may arise from, e.g., a better rejection of backgrounds). 

We also note that since [ROI, AE]o was selected using simulated data, other choices 
of [ROI, AE] may provide a stronger constraint on the data. Also, the simulated data is 
not a perfect representation of the data. Indeed, there are certain regions in the sky where 
the simulations do not model the data perfectly, and the “expected” limits using MG data 
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may differ from the limit obtained on the real data (see Appendix]^. One notable example 
is in the GC and in the Inner Galaxy region, which has led to claims of a 7 -ray excess, 
see Appendix However, an imperfect modeling of the sky does not affect the validity 
of onr constraints. We nse the simnlations as a tool to pick [ROI, AE] in an nnbiased 
way. Even if the simnlations were a totally inaccurate representation of the real data, it 
would not invalidate our limits, although other choices of [ROI, AE] would provide stronger 
constraints. 

We note that for prompt radiation we include the effects of the Fermz-LAT’s PSF, by 
performing its convolution with the J-factors, using the Fermz-LAT ScienceTools. For the 
constraints that include prompt and IGS, however, convolving the PSF for the DM signal 
calculation is computationally intensive, so we do not account for these effects. To see by 
how much this could potentially affect our limits, we constrained the ROIs to have a shape 
which is safe w.r.t. the PSF containment radius at the lowest energies considered. If the 
ROI includes a portion of GG (i.e. Ai > 0°), then we require the width of this window to be 
at least 6 ° (i.e. Ai > 3°); for the width of the top and bottom of the ROI shape (resembling 
crescents) we require that R < Ai (so the ROI is a circle), R < \/A\? -|- A^ (so the two 
crescents have no tips), and R > 4° + Ah (so the two crescents are thick enough). The 
upper bounds thus obtained are only degraded by at most ~ 20 — 40% with respect to the 
unconstrained-ROI case. This is a small number; especially in view of the fact that the 
largest uncertainty for the DM IGS signal comes from the value of the magnitude of the 
local magnetic held, see Appendix [Dj 

We note that systematic effects of the PSF are not included in our analysis, as they are 
much smaller than the other sources of systematic uncertainty considered, such as in the 
IGS signal and DM density prohle. 

3.5. Illustration of Procedure 

An illustration of our method is shown in Fig. The left plot shows the count spectrum 
from one of the MG data sets for the ROI shown in the inset. The green triangles show the 
spectrum for a 1.5 TeV DM annihilating to bb, assuming isothermally distributed DM, with 
the cross section set at the 95% G.L. upper limit. This limit is derived by requiring that 
the number of signal events in the optimal energy range from 68 GeV to 142 GeV (vertical 
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DM DM b b, /Dqm = 1.5 TeV, Isothermal 


DM DM b b, Isothermal 




Figure 3. Left: Count spectrum from one of the MC data sets for the ROI shown in the inset. 
The green points show the spectrum for 1.5 TeV DM distributed according to the Isothermal profile, 
annihilating to bb, with a cross section chosen such that the number of signal events in the energy 
range from 68 GeV to 142 GeV (vertical brown lines) is larger than the number of events in the MC 
data (at 95% C.L.), as given by Eq. (19). Since the simulated data only contains photons up to 
460 GeV, we extrapolate it to 750 GeV (red points), using a power-law fit to the photon spectrum 
above ~ 6.2 GeV. See Appendix [E] for more details. Right: The best cross-section limit averaged 
over all ten MC data sets is shown with a green solid line, while the individual cross-section limits 
for each of the 10 MC data sets are shown with dashed gray lines. As explained in 13.5, the average 
cross-section limit is used as a figure of merit for our ROI/energy range optimization. 


brown lines) be larger than the number of events in the MC data as given by Eq. (19), where 
we set a = 0.95. The number of events in this ROI and energy range will fluctuate from 
one MC data set to another, and we calculate the average cross-section limit for all ten MC 
data sets. We show the best average cross-section limit as a function of DM mass with a 
green solid line in Fig. (right), together with the cross-section limit for the ten individual 
MC data sets (dashed gray lines). In Fig. we masked all point sources and included both 
front- and back-converting events. 


We now have all the ingredients put in place for calculating constraints from the 7 - 
ray sky observed by the Fermi-LAT. In the next section we give the 95% C.L. bounds on 
the annihilation cross section (upper bound) and on the DM lifetime (lower bound) for 
annihilations and decays into various SM modes, respectively. 


16 

































4. RESULTS AND DISCUSSION 


In this section we give the results from the optimization procedure described in ^ We 
emphasize that the constraints obtained in this study are conservative and robust, since they 
do not depend on the modeling and subsequent subtraction of astrophysical foregrounds. 
In ^4.1| ( ^4.2[ ) we discuss the constraints on annihilating (decaying) DM. Additionally, in 
Appendix we use our method to derive bounds on models invoked to explain a putative 
7 -ray excess at the GC [251136] . The effect on our constraints due to different choices of 
source-masking, and due to the variation of ICS parameters is discussed in Appendix and 
Appendix respectively. 


4.1. Constraints on Dark Matter Annihilation 


The constraints on the DM-annihilation cross section as a function of DM mass are 
presented in Fig. for annihilation to and (f)cj), where (j) decays either 

only to (with rriff, = 0.1 GeV), or only to (with = 0.9 GeV), or to 

I and TT+TT” in the ratio 1:1:2 (with = 0.9 GeV). Fig. shows the results for 
the final states bb, and gg. In all cases we present the results for four different 


assumptions about the DM profile p{r) introduced in S2.1 We note that each DM mass for 
each spatial distribution and final state choice has been separately optimized, and an optimal 
ROI, ROIo,i, and photon energy range, AEo^i, were obtained to set the 95% C.L. constraint. 
In Appendix]^ we illustrate how the optimal ROI and energy range change for various DM 
density profiles and for different DM masses (see Figs. and . 


The constraints disfavor the thermal WIMP cross section for low DM masses and for 
the cuspiest profiles (mostly the NFWc profile). For those cases in which the final states 
contain high-energy electrons, i.e. Fig. there is a contribution from prompt radiation from 
FSR as well as ICS. The latter, while more uncertain, considerably strengthens the bounds, 
especially for high DM masses. In Fig. |^the shaded band denotes the constraint from ICS as 
the magnitude of the Galactic magnetic field at our Solar System’s location, Bq is varied from 
1 fiG to 10 /iG and correspondingly the diffusion coefficient Dq from 4.797 x 10^® cm^/s to 

^ Note that limits for the W'^W~ channel extend to todm < mw- In this region, the W~^W~ final state 
is not produced on-shell, but instead the annihilation is to a three- or four-body final state consisting of 
leptons and/or quarks through off-shell W^. (The expected cross-section in any concrete DM model for 
the off-shell process would be highly suppressed compared to the on-shell process.) 
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Figure 4. 95% C.L. upper limits on DM annihilation cross section vs. DM mass from Fermi-LAT’s 
inclusive photon spectrum for the indicated final states. Each plot shows constraints for the Isothermal 
(green), NEW (red), Einasto (blue), and NFWc (orange) DM density profiles. Solid lines show constraints 
from the inclusion of only the prompt radiation from the annihilation, while the bands include the ICS 
off background light, with the Galactic B-field varying within 1 — 10 yG and Dq within Do^min ~ Do,max 
(bottom-top of band). When available, we show the limits from the P7REP analysis of 15 dwarf spheroidal 
galaxies with a cyan dashed line [12]. For the XDM models we show the approximate regions (gray) in 
which annihilating DM could account for the PAMELA/Fermi/AMS-02 cosmic-ray excesses. The best-fit 
parameters from m are shown as black dots. 18 
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Figure 5. 95% C.L. upper limits on DM annihilation cross section vs. DM mass from Fermi-LAT’s inclusive 
photon spectrum for the indicated final states. Each plot shows constraints for the Isothermal (green), NFW 
(red), Einasto (blue), and NFWc (orange) DM density profiles. Solid lines show constraints derived from 
including only the prompt radiation produced in the annihilation process (i.e. final-state radiation or in the 
decay of hadrons). When available, we show the limits from the 4-year P7REP analysis of 15 nearby dwarf 
spheroidal galaxies with a cyan dashed line m- 


6.311 X lO^® cm^/s (see g2.2[). The propagation was performed as described in g2.2 


i.e. over 


a cylindrical geometry with radius Rh = 20 kpc and half-height = 4 kpc. With ICS 
included and for cuspy profiles, DM annihilation to leptonic hnal states, particularly for 
electrons, can be probed well into the annihilation-cross-section regime of a thermal relic 
that freezes out early in the Universe, (cTn)reiic ~ 3 x 10“^® cm®/s. The inclusion of extra 
particle content in DM annihilations, namely the particle (/>, is motivated by the best fit to 
the PAMELA, Fermi, and AMS-02 cosmic-ray positron and electron data [781 - I8T] . if those 
excesses are interpreted as coming from DM annihilation. Fig. shows the approximate 
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regions (shaded gray) in the cross-section-versus-niass plane, in which annihilating DM 
could offer an explanation for these excesses. These regions are meant to be illustrative only 
and chosen so that they contain the parameter choices found in [HH], shown with black dots. 
(See also |97].) The inclusion of ICS severely constrains the favored parameter regions for 
all prohles except isothermal, while including only the prompt signal challenges the favored 
regions only for the cuspy NFWc prohle. 

The constraints from |12) . which, using 4 years of P7REP data, analyzed 15 dwarf 
spheroidal satellite galaxies (dSph) of the Milky Way to set robust constraints on DM, 
are shown in Fig. |^with a cyan dashed line. Due to the dSph’s proximity, high DM content, 
and lack of astrophysical foregrounds, they are excellent targets to search for annihilating 
DM. Moreover, the available data on the velocity distribution of the stars in the dSph al¬ 
lows one to predict rather accurately the expected 7 -ray flux from DM annihilation. This 
prediction is not subject to the same uncertainties as the expected flux in the Milky Way 
halo, which suffers from large uncertainties in the DM density prohle. Our constraints are 
stronger than the dSph constraints over much of the DM mass range and for several of the 
DM prohles that we consider, especially at high energies. For DM masses < 10 GeV, our 
constraints are stronger than the dSph constraints for the NFWc prohle, and comparable in 
strength for the Einasto prohle, although weaker for the NFW and isothermal prohles. New 
results using P 8 data to perform a similar analysis of the dwarf galaxies are expected soon 
and are somewhat more stringent than the P7REP results. 

Notice that some of the ICS-inclusive limits are actually weaker than the ones with prompt 
radiation only. This might seem puzzling, as for a given ROI and energy range, the signal 
that includes prompt and ICS is obviously larger than the one with prompt only and should 
lead to more stringent constraints. However, our ROI and energy range used to derive the 
limits from the data are dictated by the optimization of the average MC limit, such that the 
optimized ROIs and energy ranges for prompt+ICS and prompt-only might diher from each 
other. If one considers this along with the fact that the simulated data sets are not perfect 
representations of the real sky, the limit that includes ICS can be weaker on occasion than 
the prompt-only limit. 

It is useful to compare our limits with those obtained from similar analyses in the liter¬ 
ature where no attempt was made to model the astrophysical foregrounds. These analyses 
usually differ in their choice of DM-pro£le parameters, their procedure for constructing the 
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limits (Gaussian error on flux versus Poisson limit on counts), their choice of propagation 
models for the ICS signal, and the data energy range utilized. Nevertheless, we can try 
to single out the effect of our ROI and energy-range optimization method alone by rescal¬ 
ing these other results to compensate for the different choices mentioned above. In iza. 
the limit was also constructed by scanning over a few differently shaped and located ROIs. 
Consequently our results are only within a factor of 1-2 stronger than theirs, across all chan¬ 
nels. In |22| . the construction of the bound is quite different from ours, and our results are 
around 2 times more stringent than theirs. In |2U) . an optimization procedure is performed 
on ROIs that look very different from ours, and a less extensive optimization is done on the 
energy window. For annihilations we improve on these limits by a factor of 1-20, depending 
on channel and profile, and by a factor 2-4 when including ICS. In [23], the ROI is opti¬ 


mized using the signal-over-background ratio as a figure of merit. For harder spectra, our 
improvement is between a factor of 3-8, while for softer spectra, the improvement is a factor 
of 1-4. 


4.2. Constraints on Dark Matter Decays 

While a favorite target for the DM annihilation rate comes from the thermal freeze-out 
of a thermal relic, which gives the correct present-day abundance, for decaying DM no 


such “favored” lifetime exists — the DM lifetime only has to be larger than the age of the 


Universe. One possible target comes from explaining the rising fraction in the cosmic-ray 
positron spectrum with DM decays to final states that produce high-energy electrons and 
positrons, with the preferred DM lifetime being in the range 10^® — 10^^ s, depending on 


the precise final states and astrophysical assumptions [211 Ea EEHinSI. Such lifetimes do 


not only have a phenomenological motivation, but also arise naturally for TeV-scale DM 
particles that decay via a dimension-six operator generated near the scale of Grand Unified 
Theories (GUT’s), M ~ 10^® GeV, namely 



(25) 


For example, in [99] DM decaying via dimension-six operators in supersymmetric GUT’s 
were posited to explain the cosmic-ray data from PAMELA. 


21 






DM -> e 


DM -> jj* jj 




DM -> T* T 


DM -> (p (p-> 2e* 2e 




DM (p (p^ 2iJ* 2iJ 


DM^(p(p^2e* 2e“, 2iJ* 2/j“, 2tt^ 27t“ 




Figure 6. 95% C.L. lower limits on DM decay lifetime vs. DM mass from Fermi-LAT’s inclusive photon 
spectrum for the indicated final states. Shown are constraints for the NFW profile (the other profiles are 
virtually identical). Solid lines show constraints derived from including only the prompt radiation produced 
in the annihilation process (i.e. final-state radiation or in the decay of hadrons), while the bands include the 
ICS off background light, with the Galactic B-field varying within 1 — 10 and Dq within Do^min ~ Do,max 
(bottom-top of band, when visible). For some models we show the approximate regions (gray) in which 
decaying DM could account for the PAMELA/Fermi/AMS-02 cosmic-ray excesses. The best-fit parameters 
from IHH] are shown as black dots. 
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Figure 7. 95% C.L. lower limits on DM decay lifetime vs. DM mass from Fermi-LAT’s inclusive photon 
spectrum for the indicated final states. Shown are constraints for the NFW profile (the other profiles are 
virtually identical). The constraints are derived from including only the prompt radiation produced in the 
annihilation process (i.e. final-state radiation or in the decay of hadrons). 


The results for DM decays to leptonic and 00 final states are included in Fig. whereas 
those decays to bb,uu, gg,W~^W~ are shown in Fig. We only show the constraints for 
the NFW profile, as the other profiles lead to virtually identical constraints. As in the case 
for DM annihilation, we include ICS for decaying DM for the leptonic final states only. 
The additional ICS component, while very sensitive to the value of the Galactic magnetic 
field, can enhance the constraints significantly, as in the case for annihilating DM. Note 
that the bounds from prompt radiation start to deteriorate near DM masses of 1.5 TeV due 
to the maximum-energy selection of 750 GeV used in this study. Our constraints compare 
favorably with existing constraints in the literature; for example, they are a factor of 2-3 
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stronger compared to [ 2 m - [ 22 ] . 


While the DM decay lifetime can span an enormous range consistent with all astrophysical 
data, there are many scenarios that are being probed by the constraints presented in this 
analysis. In particular, Fig. shows with a gray shaded parallelogram the approximate 
preferred regions in which decaying DM can explain the cosmic-ray positron and electron 
data. Black dots indicate the best-fit regions found in |98| . although note that these results do 
not include the latest data release from AMS-02 |HI] (a more careful analysis of the preferred 
regions is beyond the scope of this paper); nevertheless, we expect that the preferred regions 
would not shift significantly, and our regions are meant to be taken as a useful but rough 
qualitative guide only. We see that decays to r+r“ are thoroughly disfavored, but our 
constraints for other channels are not strong enough to probe the relevant parameter regions. 


5. CONCLUSIONS 


This paper presented a conservative method for setting constraints on 7 rays originating 
from DM annihilation and decay, which does not rely on modeling of astrophysical fore¬ 
grounds when setting a limit. Optimal regions in the sky and energy were obtained by using 
simulations of the 7 -ray sky, and a constraint was found by only requiring that the DM 
signal does not over-predict the observed photon counts. 

For models of both annihilating and decaying DM, this method allows us to constrain 
theoretically-motivated parameter regions. For example, for cuspy enough profiles (e.g., con¬ 
tracted NFW), our method is able to disfavor the thermal-relic cross section for some leptonic 
and hadronic final states. Also, for steep-enough profiles, our constraints disfavor various 
annihilating DM scenarios designed to explain the PAMELA/Fermz/AMS-02 cosmic-ray 
positron and electron data. For decaying DM, a wide range of lifetimes are excluded for 
various SM final states, including the preferred parameter regions for DM decaying to 
to explain the PAMELA/Fermz/AMS-02 data. The conservative constraints obtained in 
this study are often competitive with, and in many cases stronger than, other available 
constraints in the literature. 
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A. Constraints on DM Models invoked to explain 7 Rays from Inner Galaxy 


In this appendix we address claims made by several groups in recent years regarding a 
7-ray excess from ~ 300 MeV to ~ 5 GeV, peaking in the 1-3 GeV window, in the Inner 
Galaxy While modeling uncertainties are large and the excess may very well have 

a non-DM origin, we use our method to set constraints on DM scenarios that have been 
invoked to explain the excess. Since we perform no foreground subtraction, a priori we do 
not expect the limits derived with our method to disfavor the best-ht DM scenarios found 
in the literature; nevertheless, it is worthwhile to perform a careful check. 

The best £t for WIMP DM found in [321 133] is for ~ 30 — 40 GeV DM annihilating 
predominantly to bb. Furthermore, the spatial distribution of the putative signal is best fit 
by a generalized NFW prohle. 


pnfw, 7(^) 


_Po_ 

{r/r^y (1 + r/ry^-^ ’ 


(Al) 


with a best £t obtained for 7 1.26, although any 7 in the range ~ 1.1 — 1.4 allows for 

a reasonable £t. Analyses by other groups give results that are broadly consistent with the 
hndings in [32l [33] . In [35], it was found that DM annihilating dominantly to bh but with 
some admixture of also provides a good £t. Other annihilation channels may also be 
possible |1U6] . 


In Fig. I^we show the results of our optimization procedure applied to generalized NFW, 


Eq. (Al), with parameters chosen from best hts found in [3T| l33l l36] (which differ in part 
from the assumptions made in ^4.1 ). The authors of H (03) exclude from their analysis a 
band around the GP dehned by |&| < 1° (2°), thus not specifying a specihc DM distribution 
within this latitude. We therefore use our usual ROIs shown in Fig. but mask a square 
centered on the GG of side 2° (4°). We show DM annihilating to bb (left plot) and 
(right plot). Unsurprisingly, the bounds that we obtain on the annihilation cross section 
are still a factor of ~ 3 or more from probing the best-ht regions shown with open or closed 
contours in Fig. As a reference for the reader, adopting all the assumptions in [33], for 
annihilation into bh, and choosing rriDM = 25 GeV, the optimal ROI found with our method 
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Figure 8. 95% C.L. annihilation cross section upper limits on DM annihilating to bb (left) and 
T^T~ (right) for an NFWc profile with various inner slopes and local DM densities (note that the 


assumptions made in deriving these limits differ in part from those made in ^.1). Also shown are 


the preferred regions from [311 [331136] for DM to fit the claimed Galactic-Center 7 -ray “excess”. The 
constraints have been computed with the same model assumptions as the best-fit regions (including 
masking a square centered on the GC of side 2 ° or 4° for analyses that excluded a band around the 
GP with the same thickness - see text for details). We also show with a cyan dashed line the limit 
obtained from the 4-year P7REP analysis of 15 nearby dwarf spheroidal galaxies |12j . 


is determined by the following parameters: R = 2°, Ab = 1.98°, Ai = 0.54°, while the 
optimal energy range is 1.9 GeV< E < 4.0 GeV. 


B. Dependence of Optimal ROI and Energy Range on DM Profile and DM Mass 


The optimal ROI and photon-energy range are found separately for each choice of DM 
spatial distribution, mass, and final state. In this section, we briefly illustrate the generic 
features of the optimal search region and its dependence on the theory hypothesis. Fig. 
shows the obtained ROI and energy range for DM annihilation to bb for each of the four 
spatial distributions studied, and for a fixed DM mass of 25 GeV. For this final state, with 
the exception of NFWc, where it is beneficial to look near the GG, the optimal regions 
in the sky involve semi-circular regions, symmetric in latitude b, with the GG removed. 
Furthermore, we find narrower optimal energy ranges for NFWc-distributed DM. 


For the bb final state, the effect of varying the DM mass is addressed in Fig. [1^ where the 
optimal regions are shown for two different masses: 350 GeV and 7 TeV, assuming NFWc- 
distributed DM. As the DM mass is increased, the strongest optimal regions are obtained 
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Figure 9. Count spectrum for 25 GeV DM annihilating to bb for various DM density profiles. The vertical 
(brown) lines show the optimal energy range for each DM model assumption. The inset shows the optimal 
ROI. Note that PSF-convolution effects were included for the DM signal. The quoted {av) is the annihilation 
cross section that saturates the 95% C.L. from the data. 


by including semi-circular regions in latitude, in addition to a rectangular area around the 


GC. We note that finite-resolution effects were included, by convolving the instrument’s PSF 


with the J-factors, in the DM signal for all of the results in 


Fig. [^and Fig. [TO 


C. Effect of Source Masking and Choice of Pront-/Back-converting events on Limits 


In this appendix we investigate the effect on the DM-cross-sect ion upper limits when 
masking known point sources and using front- and/or back-converting events. 

Masking known sources reduces the observed counts in an ROI and can strengthen the 
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DM DM b b, /Ddm = 7 TeV, Einasto 


DM DM bb, /Ddm = 350 GeV, Einasto 



Figure 10. Count spectrum for 350 GeV (left) and 7 TeV (right) DM annihilating to bb, assuming an 
Einasto profile. The vertical (brown) lines show the optimal energy range for each DM-model assumption. 
The inset shows the optimal ROI. Note that PSF-convolution effects were included for the DM signal. The 
quoted (av) is the annihilation cross section that saturates the 95% C.L. from the data. 

DM constraints, assnming that the masking does not also remove mnch of a potential DM 
signal. This is the case if the ROI is large, as it is expected to be for decaying DM, or for 
annihilating DM with shallow DM density prohles (e.g., isothermal). Since astrophysical 
point sonrces at very large energies (> 20 GeV) typically exhibit a small flnx, their masking 
is expected to improve the limits for lower DM masses. For very cnspy prohles the ROIs tend 
to be small and concentrated around the GC region, where the number of known sources is 
also large; in this case, masking all the point sources would remove most of the DM signal 
as well and will thus not likely lead to stronger limits. 

The amount of sky that needs to be masked to remove a point source depends on the 
Fermi-LAT PSF, which depends on the energy and on where the photon converts in the 
detector. In particular, photons that convert to an e’''e“ pair in the front part of the Fermi- 
LAT (consisting of the hrst 12 layers of thin tungsten foil) have a better angular resolution 
(smaller PSF) than those photons that convert in the back (next 4 layers of thick tungsten 
foils). For very cnspy prohles the choice of including only front- or only back-converting 
events, or both, could potentially have important effects on the constraints. 

We obtain the point-source coordinates from the 3FGL catalog |95] and exclude all the 
photons contained in pixels whose center lies within an angular radius of 2 6q8{E) from any 
point source; here 6q8 is an approximation of the energy-dependent P7REP 68% point-source 
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Figure 11. Ratio of expected cross section upper limits vs. DM mass from simulated MC data 
for DM annihilation to bb for isothermal (left) and NFWc (right) profiles. The denominator of 
the ratio, ((Tn)F+B,s.m.) is the cross section upper limit obtained when masking all known point 
sources in the 5-year Fermi-hAT point-source catalog, outside a 2° x 2° square centered at the GC 
and including front- and back-converting events. The numerators of the ratios, {crv)i, are the cross 
section upper limits obtained when masking all known point sources outside the 2° x 2° GC square 
(solid lines), masking only those sources detected at more than lOcr (outside the same 2° x 2° GC 
square) (dashed lines), and masking no sources (dotted lines). In each case we either include both 
front- and back-converting events (blue lines), only front-converting events (red lines), and only 
back-converting events (green lines). 


containment angle, 


968(£)ri = \/c2(B/lGeV)-"'’ + c5, 


(Cl) 


and the parameters for (front-, back-) converting events are cq = (0.645,1.103), ci = 
(0.0821,0.166), and f3 = (0.762,0.750). No source masking is performed within the in¬ 
ner 2° X 2° square at the GC, where the density of sources is very high and the expected 
DM signal peaks. 

The effect on the cross-section upper limits versus DM mass, when masking known point 
sources, and when including front- and/or back-converting events, is shown in Fig. 11 on 


simulated data sets. The left (right) plot assumes DM annihilation to bh for our choice of 
an isothermal (NFW-contracted) density prohle. We choose a shallow and cuspy prohle to 
see how the results depend on having either large or small optimized ROIs, respectively. 
For each DM mass, and for each choice of source masking and inclusion of front-/back- 
converting events, we optimize the ROI choice and derive the average limit obtained from 
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the ten simulated MC data sets. In Fig. [m we show a ratio of expected cross section upper 
limits versus DM mass: the denominator of the ratios, (crn)F+B,s.m .5 is the cross section upper 
limit obtained when masking all known point sources as described above and including front- 
and back-converting events; the numerators of the ratios, (crn)j, are the cross section upper 
limits obtained when masking all known point sources (outside the 2° x 2° GC square) (solid 
lines), masking only those sources detected at more than lOa (outside the same 2° x 2° GC 
square) (dashed lines), and masking no sources (dotted lines). In each case we either include 
both front- and back-converting events (blue lines), only front-converting events (red lines), 
or only back-converting events (green lines). 


We see from Fig. IG that, at least for the two annihilation models considered in this 
section, the expected limits are the same within 0(10 — 30%). Moreover, the strongest 
constraints are generically obtained when masking all point sources. For DM masses below 
~ 50 GeV and cuspy prohles, the inclusion of only front-converting events is expected to 
provide the strongest constraints, but only marginally so. Above ~ 50 GeV, the inclusion of 
both front- and back-converting events is best, since the photons produced in the annihilation 
of DM have such high energies that the PSF effects are negligible, and the inclusion of as 
much data as possible leads to stronger expected limits. 

Based on this, we conclude that the effect of source masking and choice of front-/back- 
converting events is not large on our results. We also note that the inclusion of both event- 
conversion types and the masking of point sources (blue solid line in Fig. [II] ) is expected 
to give constraints that are among the best. We thus make this our standard choice when 
showing the results in ® 


D. Inverse Compton Scattering 

In this appendix we discuss how the results from ^ depend on the parameters in the 
IGS computation performed in GALPROP. The amount of IGS radiation depends sensitively 
on various key propagation parameters whose values are not known to a satisfactory degree. 
Here we describe the effect on our constraints from varying these parameters in order to 
capture some of the systematic uncertainties associated with the DM-generated IGS signal. 

We study how different models of propagation impact our results. We use, as a start¬ 
ing point, the Fermz-LAT results from IMl, in which various propagation models are £t 
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Figure 12. Ratios of cross-section upper limits from simulated data on DM annihilation to e'’'e“ 
for an NFW profile, including prompt and ICS radiations, for different values of the Galactic 
magnetic field (left) and different combinations of other propagation model parameters (right). 
The magnetic field has the largest effect on our analysis. 


to cosmic-ray spectra for various choices of the region of containment of the cosmic rays 
(parametrized with a cylindrical geometry of half-height Zh and radius Rh)- In our study, 
we vary Zh and Rh, and two other important parameters that have a big effect on the DM 
ICS signal, namely the Galactic magnetic field value in the Solar System, Bo, and the spatial 
diffusion coefficient Do- The values used in our study are: 

1. Zh = 4, 6, 8,10 kpc 

2. Rh = 20, 30 kpc 

3. Dq = Do,min, Dg^max, where Do,min and Do,max are the minimum and maximum values 
of Do spanned by the various GALPROP models studied in |107] for a given {zh, Rh)- 


4. Do = 1,5,10 fiG 


As an illustration of the dependence of the DM ICS signal on these parameters, Fig. 12 


shows the constraint on DM annihilation to assuming an NFWc DM profile. The 

greatest effect on the uncertainty of the DM ICS signal originates from the variation in the 
magnitude of Dq, as clearly shown in the left plot. Varying the other parameters (right plot) 
has less of an effect on the DM ICS signal. We are therefore allowed, when showing the 
results in ® to fix Zh = 4 kpc and Rh = 20 kpc; whereas we show the variation of our results 
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with Bq and correspondingly Dq = -Do,min = 4.797 x 10^® cni^/s for Dq = 1 fiG (parameters 
yielding the strongest constraints) and Dq = Do,max = 6.311 x 10^® cm^/s for Bq = 10 /iG 
(parameters yielding the weakest constraints). 


E. Details on the Simulated data sets 


The optimization procednre described in 53.4 to find the optimal ROIs and energy ranges, 


[ROI, AE]o, is performed on ten simnlated data sets, each a 5.84-year representation of the 
7 -ray sky. Here we provide a few more details on the simnlations. 

The generation of mock Fermi-LAT observations was carried ont with the gtobssim 
rontine, part of the Fermi Science Tools package v9r34pl. Its ontpnt is a list of MC- 
simnlated y-ray events with relative spatial direction, arrival time and energy, distribnted 
according to an inpnt sonrce model and IRFs. 

A nnmber of model elements were pnt into gtobssim (see |in8) ). These inclnde the Fermi- 
LAT Collaboration’s model of the diffnse Galactic component]^ the isotropic component 
(derived for Pass 7 Reprocessed Clean front and back IRFs) and the 3FGL sonrce catalog 
for point and small extended sonrces |95) . 

In addition, the fnll-sky simnlations were calcnlated throngh gtobssim with the actnal 
pointing and livetime history (FT2 file) of the Fermi-LAT for the first 5.84 years of the 
scientific phase of the mission. The sonrce model simnlated did not contain the Earth’s 
Limb emission, which is negligible at energies above 1 GeV, compared to the celestial y-ray 
signal, when a zenith angle < 100° cnt is applied. The gtobssim tool convolves the flnx 
components mentioned with the Fermz-LAT’s response, i.e. PSF, energy dispersion, and 
effective area. 

Ten instances of the MG gtobssim-generated data were rnn, each with an independent 
starting seed and the same sonrce model; thns obtaining ten statistically independent in¬ 
stances of the y-ray sky. The same event selection criteria were nsed for the MG data sets 
as for the real data. One important difference between the simnlated data sets and the real 
data is the energy range. Each simnlated data set was calcnlated in an energy range of 
0.5 GeV to 500 GeV (as opposed to 1.5 GeV to 750 GeV for the actnal data). The npper 
bonnd of 500 GeV in the gtobssim simnlations is the npper limit in the energy map of the 


gll_iem_v05_revl.fit 

® iso clean front v05.txt and iso clean back v05.txt 
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Figure 13. Comparison between average MC-based expected (dashed) and real-data (solid) 95% 
C.L. annihilation cross section npper limits on DM annihilating to bb (left) and e'’'e“ (right) for the 
Isothermal (green), NFW (red), Einasto (bine), and NFWc (orange) DM density profiles (we only 
consider prompt photons). The popnlation standard deviations of the limits from the 10 individnal 
MC simnlations are also shown as shadings aronnd the dashed lines. 


interstellar diffuse model Hffll. To deal with this mismatch, we simply ht a power-law curve 
to each of the ten simulated data spectra for 6.2 GeV< E < 460 GeV that we obtain for 
each ROI, and extrapolate it to 750 GeV. (The lower value of 6.2 GeV is low enough to 
have enough photons to perform a meaningful ht even for small ROIs, and high enough for 
a single featureless power law to provide a reasonable ht to the spectra. The upper value of 
460 GeV is low enough to avoid count leakages due to hnite energy resolution on the sharp 
500 GeV input-energy cutoh.) We then populate each bin above 460 GeV with a random 
number of events chosen from a Poisson distribution whose expectation value equals the 
extrapolated value in a given bin. The subsequent optimal ROI and energy range for each 
theory hypothesis Th is found using the original plus extrapolated spectra. 


F. Comparison of limits between simulated and real data 


In this appendix we compare the results derived from the real data with those derived 
from simulated data. Since our simulated data is of course not a perfect representation of 
the real data, we do not expect that the limits derived on the real data will agree perfectly 
with the limits derived on simulated data. 


Fig. 13 compares the simulated and observed limit on DM annihilation to bb (left plot) 
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and e+e“, including only prompt photons, (right plot), for the four different DM density 
profiles introduced in § 2 . 1 [ Since the simulated data used in this study consists of 10 
statistically independent realizations of the 7 -ray sky, we present the arithmetic mean of 
the 10 limits (dashed lines) and the standard deviation of the population (shaded bands), 
as well as the observed limits (solid lines). We see that the limits derived using real versus 
simulated data agree over a wide range of masses and profiles. 
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